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Abstract 

Recent theoretical results establish that time- consistent valuations (i.e. pricing operators) 
can be created by backward iteration of one-period valuations. In this paper we investigate 
the continuous-time limits of well-known actuarial premium principles when such back- 
ward iteration procedures are applied. We show that the one-period variance premium 
principle converges to the non-linear exponential indifference valuation. Furthermore, we 
study the convergence of the one-period standard-deviation principle and establish that the 
Cost-of-Capital principle, which is widely used by the insurance industry, converges to the 
same limit as the standard- deviation principle. Finally, we study the connections between 
our time -consistent pricing operators, Good Deal Bound pricing and pricing under model 
ambiguity. 
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1 Introduction 



Standard actuarial premium principles usually consider a static premium calculation problem: 
what is the price today of an insurance contract with payoff at time T. See, for example, the text- 
books by Buhlmann (1970), Gerber (1979), or Kaas et al. (2008). Also, the study of risk measures, 
and the closely related concept of monetary risk measures has been studied in such a static 
setting. See, for example, Artzner et al. (1999), Cheridito et al. (2005). Also the study of utility in- 
difference valuations has mainly confined itself to this static setting. For different applications 
we mention a few papers: Henderson (2002), Young and Zariphopoulou (2002), Hobson (2004), 
Musiela and Zariphopoulou (2004), Monoyios (2006), and the recent book by Carmona (2009). 

Financial pricing usually considers a "dynamic" pricing problem: how does the price evolve 
over time until the final payoff date T. This dynamic perspective is driven by the focus on hedg- 
ing and replication. This literature started by the seminal paper of Black and Scholes (1973) 
and has been immensely generalised to broad classes of securities and stochastic processes, see 
Delbaen and Schachermayer (1994). 

In recent years, researchers have begun to investigate risk measures in a dynamic setting, 
where the question of constructing time-consistent (or "dynamic") risk-measures has been in- 
vestigated. See, Riedel (2004), Cheridito et al. (2006), Roorda et al. (2005), Rosazza Gianin (2006), 
Artzner et al. (2007). In a recent paper lobert and Rogers (2008) show how time-consistent valu- 
ations can be constructed via backward induction of static one-period risk-measures (or "valu- 
ations"). 

In this paper we want to investigate well-known actuarial premium principles such as the 
one-period variance principle and the standard- deviation principle, and study their time-consistent 
extension. The method we use to construct these extensions is to first consider one-period 
valuations, then extend this to a multi-period setting using the backward iteration method of 
Jobert and Rogers (2008) for a given discrete time-step At, and finally consider the continuous- 
time limit for A t — ► 0. We show that the one-period variance premium principle converges to the 
non-linear exponential indifference valuation. Furthermore, we study the convergence of the 
one-period standard- deviation principle and establish that the Cost-of-Capital principle, which 
is widely used by the insurance industry, converges to the same limit as the standard-deviation 
principle. 

The rest of this paper is organised as follows. We start in Section 2 where we focus initially on 
the case of pure insurance risk. This allows us to clearly demonstrate the construction we use 
to derive the continuous-time limit of time-consistent actuarial pricing principles. We derive 
the limit for the Variance principle, the Standard-Deviation principle and the Cost-of-Capital 
principle. . . . We summarise and conclude in Section 6. 
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2 Variance Pricing 



2.1 Diffusion Setting 

We start by considering an unhedgeable insurance process y t , which is given by a diffusion equa- 
tion: 

dy = a{t,y)dt + b{t,y)dW. (2.1) 

We also consider a discretisation scheme for the insurance process in the form of a binomial 
tree: 

f +b\TKi withprob. I 
y(t+At) = y(t) + aAt + \ r— v \ (2.2) 

{ -bVAt withprob.^ 

where we have suppressed the dependence of a and b on (f,y) to lighten the notation. Note 
that we have restricted ourselves to a Markovian diffusion setting, which allows us to give a 
very simple mathematical derivation of our results. In Section 4 we discuss alternatives for this 
restrictive assumption. 

Given the discrete-time setting (2.2), we can now create time-consistent pricing operators 
("valuations"), using the backward induction method of Jobert and Rogers (2008). Let us denote 
a one-step valuation by n [] , and the resulting price at (t,y) by n(t,y): 

n(t,y) = U t [n(t + At,y{t + At))] (2.3) 

In words: the price n(t,y) is obtained by applying at time t the one-step valuation U t [] to the 
random variable n[t + At, y{t + At)), which is the price obtained in the previous time-step. 



2.2 Variance Principle 

If we consider an insurance contract with a payoff at time T defined as a function f{y{T)), then 
the actuarial Variance Principle Yl v t [] is defined as (see, e.g. Kaas et al., 2008) 

U v t [f(y(T))] = E t [f(y(T))] + |aVar t [/(y(T))], (2.4) 

where E f [] and Var f [] denote the expectation and variance operators conditional on the infor- 
mation available at time t under the "real-world" probability measure P. 

Note that in the standard actuarial literature (see, e.g. Kaas et al., 2008), discounting is usu- 
ally ignored. To facilitate the discussion, we will first derive the continuous-time limit of the 
variance principle without using discounting in Section 2.2.1. We will then consider case with 
discounting in Section 2.2.2, and discuss further generalisations in Section 4. 



2.2.1 No Discounting 

In the binomial tree discretisation we can obtain an explicit expression for a one-step variance 
price 7r v (f,y) by substituting (2.4) into (2.3): 

7T v [t,y(t)) = E t [n v [t + At,y{t + At))] + ±aY<m t [7T v {t + At,y(t + At))]. (2.5) 
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We are now interested in considering the limit for At — ► 0. We assume that n v {t + At, y) is suffi- 
ciently smooth to be twice continuously differentiable in y, such that we can apply for all values 
of y the Taylor approximation 

n v {t + At,y{t) + h) = n v (t + At,y{t)) + 

7i v y (t + At,y(t))h+ \n yy {t + At,y{t))h 2 + 0(h 3 ), (2.6) 

where subscripts on tt v denote partial derivatives. If we substitute this Taylor approximation for 
the binomial approximation (2.2) into (2.5) and gather all terms in ascending orders of At, we 
obtain 



n v [t,y(t))-n v [t + At,y[t)) = an y [t + At,y{t))At + 

\b 2 n yy {t + At,y{t))At+\a{bn y {t + At,y{t)]^ At + &{At 2 ). (2.7) 

If we divide by A t we obtain 

7T v (t,y{t))-n v [t + At,y{t)) , 
— ^ — - - — - = an y { t + A t, y ( t)) + 

jb 2 n yy {t + At,y{t)) + \a[bTi y {t + At,y{t))f +©{At). (2.8) 

We can now take the limit for At — ► 0. The left-hand side of (2.8) converges to (minus) the partial 
derivative of n v Q with respect to t (i.e. -n v t {t, y(f))), and we obtain 

n v t + aTi y + \b 2 n v yy + \a{bn y ) 2 = Q, (2.9) 

where we have suppressed (again) the dependence on t and y to lighten the notation. 

Note, that equation (2.9) is a semi-linear partial differential equation (pde) that describes the 
behaviour of the variance price n v {t, y) as a function of t and y. The pde is subject to the bound- 
ary condition n v (T, y[T)) = f(y{T)) which is the payoff of the insurance contract at time T. The 
existence of solutions of semi-linear pde's has been studied in the context of Backward Stochas- 
tic Differential Equations (BSDE's), and we discuss this subject further in Section 4. 

In this particular case, we can construct the solution of (2.9) explicitly by employing a Hopf- 
Cole transformation of the solution that removes the non-linearity from the pde. Consider the 
auxiliary function h v (t, y) := exp{an v (t, y)}. The original function n v (t, y) can be obtained from 
the inverse relation n v {t, y) = ^ In h v (t, y). If we now apply the chain-rule of differentiation, we 
can express the partial derivatives of n v {) in terms of h v Q as 

v iK ih y i ^ v -^> 2 f2im 

a W y a h v yy a (h v y 

If we substitute these expressions into (2.9), the non-linear terms cancel and we obtain a linear 
pde for h v {t,y): 

h v t + ah y + lb 2 h yy = 0. (2.11) 
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Hence, by considering the transformed function h v {t, y) we have managed to obtain a linear pde 
for h v (). The boundary condition at r is given by h v {T, y{T)) = exp{an v (T,y(T))} = exp{af(y(T))}. 
Using the Feynman-Kac formula, we can express the solution of (2.11) as 



h\t,y) = -E t 



,af(ym) 



y(t) = y 



(2.12) 



where the expectation is taken with respect to the stochastic process y{t) defined in equation (2.1) 
conditional on the information that at time t the process y(t) is equal to y. From the represen- 
tation (2.12) follows immediately that we can express n v {t, y) as 



7r v U,y) = -lnE f 
a 



,af{y[T)) 



yit) = y 



(2.13) 



Note that this representation of the variance-price n v () is identical to the exponential indiffer- 
ence price which has been extensively studied in recent years. See, for example Henderson 
(2002), Young and Zariphopoulou (2002) or Musiela and Zariphopoulou (2004). For an overview 
of recent advances in indifference pricing, we refer to the book by Carmona (2009). 

To summarise this section, we have established that the continuous-time limit of the iterated 
actuarial variance principle is the exponential indifference price. 



2.2.2 With Discounting 

Up to now we have ignored discounting in our derivation. (Or equivalently, we assumed that 
the interest rate is equal to zero.) In a time-consistent setting, it is important to take discounting 
into consideration, as money today cannot be compared to money tomorrow. 

If we consider the definition of the variance principle given in (2.4), it seems that we are 
adding apples and oranges together. The first term E f [/(y(T))] is a quantity in monetary units 
(say €) at time T. However, the second term Var f [f{y(T))] is basically the expectation of f{y{T)) 2 , 
and is therefore a quantity in units of (€) 2 . The way to rectify this situation is by understanding 
that the parameter a is not a dimensionless quantity, but is a quantity expressed in units of l/€. 
This should not come as a surprise. The parameter a is in fact the absolute risk aversion pa- 
rameter introduced by seminal paper by Pratt (1964) where he derives the variance principle as 
an approximation "in the small" of the price that an economic agent facing a decision under 
uncertainty should ask. 

To stress in our notation the units in which the absolute risk aversion a is expressed, we 
will rewrite the absolute risk aversion as the relative risk aversion y (also introduced by Pratt, 
1964), which is a dimensionless quantity divided by a benchmark wealth-level X(T), which is 
expressed in € at time T. If we now assume a constant rate of interest r, we can then set our 
benchmark wealth as X(T) = X e rT . Hence, we rewrite our variance principle as 

iW(ym)] = E f [/(y(T))] + l—L—Vai t [f{y{T))]. (2.14) 

X e rl 

Note, that U][] leads to a "forward" price expressed in units of € at time T. 
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Given the enhanced definition (2.14) of the variance principle including discounting, we can 
now proceed as in Section 2.2.1. For a single binomial step, we obtain the following expression 
for the price: 



n v (t,y(t)) = e-'' At \E t [n v (t + At,y(t + At))] + 



b P nt + ^ Y ^ v ( t+At 'y( t+At M- (2 - 15) 



X e 



Note that we have included an additional discounting term e~ rAf to discount the values from 
time t + A t back to time t. Using a similar derivation as before, we arrive at the following partial 
differential equation for 7r v (f, y): 

TiZ+anZ + Wn\ v + \—^—{bnl) 2 - m v = 0. (2.16) 

y z yy ^ X^e' 1 y 

This non-linear pde can again be linearised by considering the transformation h v (t,y) = exp{ y r rt n v (t, y)}, 
which leads to the following expression for the solution of (2.16): 



r 



r 



y(t) = y 



(2.17) 



From this result we see that the discounting is incorporated into the non-linear pricing formula, 
by expressing all units relative to the "benchmark wealth" X(t) = X e rt } See chapter by MuZa 
in Carmona book?? 

2.2.3 Current price as benchmark 

In the previous subsection we have taken the benchmark wealth to be a risk-free investment 
X e rt . Another interesting example is when we consider the current price n{t,y) as the bench- 
mark wealth. 

This then leads to a new pricing operator, which we will denote by n p Q. The one-step valua- 
tion is then given by 

r\t( ^ i Var f [7r p ^ + Af, y(t + Af))h 

1 y ' \ ' JK JJ> 2 ' E t [7iP{t + At,y(t + At))] I 

Hence, we assume that we want to measure the variance of n p () relative to the expected value of 
7T P (). Obviously, this is only well-defined if 7r p (f, y) is strictly positive for all (f,y). 

If we employ our Taylor-expansion and take the limit for A t — ► we obtain the following pde 

n p t + an p y + \b 2 n p y + \^{bn p f - rn p = 0. (2.19) 



J For general results concerning "benchmark pricing" in a linear setting we refer to Platen (2006) and the book 
by Platen and Heath (2006). 
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Again, we can study the solution of (2.19) by employing a transformation of the solution that 
removes the non-linearity from the pde. Consider the auxiliary function h p {t,y) := {n p {t,y)) llq . 
The original function can be obtained from the inverse relation 7i p {t, y) = {h p {t, y)) q . If we now 
apply the chain-rule of differentiation, we can express the partial derivatives of n p {) in terms of 
h p {) as 

Ti* = qthP) q - l h v , n p = q{h p ) q - l h p , n yy = qih^ 1 [^-itff + rfy) ■ (2-20) 
If we substitute these expressions into (2.19) and simplify, we obtain 

hP -(^) 2 + ^ y )-^ P = 0. (2.21) 



y 



If we choose q = 1/(1 +y), then the non-linear terms cancel and we obtain a linear pde for 

hP{t,y): 

h p t +ah p y + lb 2 h p yy -ra + T)h p = 0. (2.22) 

The boundary condition at T is given by h p (T,y{T)) = n v (T,y(T)) 1+r = f[y[T)) 1+r . Using the 
Feynman-Kac formula, we can express the solution of (2.22) as 

hP(t,y) = E t [e- ni+ ? ]iT - t] f(y(T)) l+ ?\y(t) = y], (2.23) 

where the expectation is taken with respect to the stochastic process y(f) defined in equation (2.1) 
conditional on the information that at time t the process y(f) is equal to y. From the represen- 
tation (2.23) follows immediately that we can express n p (t,y) as 

n p {t,y) = e~ r[T - t] (E f [f{y{T)) l+ ^\y{t) = y})^ , (2.24) 

Note that this representation of the price n p () arises also in the study of indifference pricing 
under power-utility functions, and the related notion of pricing under so-called "^-optimal" 
measures. See, for example Hobson (2004) and Henderson and Hobson (2009). 



2.3 Mean Value Principle 



The examples we gave in the previous subsections, are all special cases of the Mean Value Prin- 
ciple, which is defined as 



n™[f(y{T))] = v-'lEAvifiyiTM) 



(2.25) 



for any function vQ which is a convex and increasing (see Kaas et al., 2008, Chap. 5). 

Once more, we have to pay attention to units. If we want to apply a general function vQ to a 
value (expressed in units of €), we have to make sure that the argument of vQ is dimensionless. 
The easiest way to achieve this, is to express the argument for v () in "forward terms". For a single 
binomial step, we therefore obtain the following expression for the price: 



n m [t,y(t)) 



,rt 





i 




V 


V 





n m {t + At,y{t + At)) 

e r(t+At) 



(2.26) 
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We can rewrite this definition as 



n m {t,y(t)) 



,rt 



7T m [t + At,y{t + At)) 

e r(t+At) 



(2.27) 



and from this expression it is immediately obvious that the "distorted" value v[n m (t,y)/ e r ) is 
linear and therefore satisfies the Feynman-Kac formula. Which corresponds exactly to the solu- 
tions we found in the previous subsections. 

In this case, we want to go in the opposite direction and we seek the corresponding pde for 
the price n m (t,y). To do this, we will use our Taylor-expansion derivation, but with an additional 
twist: we must expand the functions vQ and as well. In particular, we seek to expand the 
function v~ l {v{x) + h) for small h. Using the identity v~ l {v{x)) = x we obtain 



v l {v{x) + h) = x + 



h 



v"(x) 
v'{x) ~ 2{v'(x)) 3 



h 2 + 0(h 3 ). 



(2.28) 



Combining this result with the Taylor expansions for vQ and n m Q on a single binomial time-step, 
and taking the limit for At -> leads to the pde: 



7T t +CITI,, +;fl 71,,,,+ ~ 7-[D7l v ) = I), 

f V 2 yy 2 y/^mfj V 



(2.29) 



where n mf {t,y) := n m {t, y)le rt is the price expressed in forward terms. 

We see in equation (2.29), that the coefficient in front of the non-linear term can be identified 
as the "local risk aversion" induced by the function v Q at the current value 7r mf () . Note, that since 
the function v() is increasing and convex by assumption, v"Q I v'Q is positive. 



3 Standard-Deviation Pricing 

3.1 Standard-Deviation Principle 

Another well-known actuarial pricing principle is the Standard Deviation Principle, (see Kaas et al., 
2008) defined as 

n s f [f(y(T))] = Er[/(y(T))] + Py/Vai t [f{y{T))]. (3.1) 

Please note that also in this case we have to be careful about the dimensionality of the parameter 
(5. Even though both the expectation and the standard deviation are expressed in units of €, the 
standard deviation and the expectation have different "time-scales". If we go down to small 
time-scales (as we will be doing when considering the limit for A t — ► 0) then due to the diffusion 
term d W of the process y, we have the property that the expectation of any function /(y) scales 
linearly with A t, but the standard deviation scales with \/a7. This means that for small A t the 
standard deviation term will completely overpower (literally!) the expectation term. Therefore, 
the only way to obtain a well- defined limit for At — ► is if we take fis/At as the parameter for the 
standard deviation principle for a binomial step. 
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Another way of understanding this result, is to consider the following example. If we want 
to compare a standard deviation measured over an annual time-step with a standard deviation 
measured over a monthly time-step, we have to scale the annual outcome with \/l/12 to get a 
fair comparison. 

Given the above discussion on the time scales, we get for a single binomial step, the following 
expression for the price: 

n%t,y(t)) = e~ r ^^ t [n%t + At,y(t + At))] + /5v^yjvav t [n s {t + At,y(t + At))]^. (3.2) 

Using a similar derivation as in Section 2.2, we arrive at the following partial differential equation 
for n s {t,y): 

n s t + an s y + \b 2 n s yy + 0yj {bn s y ) 2 - rn s = 0, (3.3) 

which can be rewritten as 

n\ + an y + \b 2 n yy + pb\7T y \ - rn s = 0. (3.4) 

This is again once more a non-linear pde. However, the non-linearity is much more benign in 
this case. Whenever the partial derivative n y does not change sign on the whole domain of y 
(i.e. the function n s is either monotonically increasing or monotonically decreasing in y), then 
(3.4) reduces to the linear pde: 

n s t + (a + pb)n s y + \b 2 n s yy -r n s = 0, (3.5) 

where the sign of ±fib depends on the (uniquely defined) sign of n y . 

Using the Feynman-Kac formula, we can represent the solution to (3.5) as: 

n\t,y) = E? [f{y{T))\y{t) = y] , (3.6) 

where E ^ [] denotes the expectation at time t with respect to the "risk- adjusted" process y de- 
fined as 

dy= [a{t,y)±pb(t,y))dt+b{t,y)dW s . (3.7) 

The drift-rate is adjusted upwards (a+fib) if the payoff /(y) is monotonically increasing in y, and 
adjusted downwards (a - fib) if /(y) is monotonically decreasing in y. So, the risk- adjustment is 
always in the "upwind" direction of the risk, thus making the price n s more expensive than the 
real-world expectation E[/(y)]. 

3.2 Cost-of-Capital Principle 

Another actuarial pricing principle is the Cost-of-Capital Principle. This was introduced by the 
Swiss insurance supervisor as a part of the method to calculate solvency capitals for insurance 
companies (Keller and Luder, 2004). 2 The Cost-of-Capital method has been widely adopted by 

2 For a critical discussion on the risk-measure implied by the Swiss Solvency Test we refer to 
Filipovic and Vogelpoth (2008) . 
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the insurance industry in Europe, and has also been prescribed as the standard method by the 
European Insurance and Pensions Supervisor for the Quantitative Impact Studies (see EIOPA, 
2010). 

The Cost-of-Capital is based on the following economic reasoning. We first consider the 
"expected loss" E[/(y(T)] of the insurance claim f[y(T)) as a basis for pricing. But this is not 
enough, the insurance company also has to hold a capital buffer against the "unexpected loss". 
This buffer is calculated as a Value-at-Risk over a time-horizon (typically 1 year) and a prob- 
ability threshold q (usually 0.995, or even higher). The unexpected loss is then calculated as 
VaRq [/(y(T)) - E[/(y(T))] ] . The capital buffer is borrowed from the shareholders of the insur- 
ance company (i.e. the buffer is subtracted from the surplus in the balance sheet). Given the 
very high confidence level, in most cases the buffer can be returned to the shareholders, how- 
ever there is a small probability (1 - q) that the capital buffer is needed to cover an unexpected 
loss. Hence, the shareholders require a compensation for this risk in the form of a "cost-of- 
capital". This cost-of-capital needs to be included in the pricing of the insurance contract. If we 
denote the cost-of-capital by 8, then the Cost-of-Capital Principle is given by 



n c t [f(ym)] = ^ t [f(ym)]+8VaR q , t f(y(T)) - E f [f(y(T))] 



(3.8) 



Please note that also in this case we have to be careful about the dimensionality of the different 
terms. First, we are comparing Value-at-Risk quantities at different time-scales, and these have 
to be scaled back to a per annum basis, to do this we divide the VaR-term by \fXi. Then, we 
must realise that the cost-of-capital 8 behaves like an interest rate: it is the compensation the in- 
surance company needs to pay to its shareholders for borrowing the buffer capital over a certain 
period. The cost-of-capital is expressed as a percentage per annum, hence over a time-step At 
the insurance company has to pay a compensation of 8 A t per € of buffer capital. As a result, we 
obtain a "net scaling" of 8 A tl y/At = <5\/a7. Note, that this is the same scaling as for the standard 
deviation principle. 

For a single time-step, we therefore get the following expression for the cost-of-capital price: 



n c (t,y(t)) = e- rAt \E t [n c (t + At,y(t + At))) + 



8VAtVaR qit n c {t + Af,y(r + Af)) - E f [n c {t + At,y{t + At))] 



(3.9) 



In the previous sections we used a binomial discretisation of the process y. However, in this 
case we have to be a bit more careful. Since we are considering a (1 - g)-quantile with very small 
probability (1 - q), a simple binomial tree approximation is too crude to obtain an accurate 
representation of the (1 - g)-quantile of the process y. We therefore consider a "quadrinomial" 
tree, where we make sure that we match the mean, the variance and the (1 - g)-quantile of 
the process y over a. At time-step. We can do this, for example, by considering the following 
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discretisation: 



y(t + At) = y{t) + aAt+ < 



+ kb 
+ lb 
-lb* 
-kb 



At with prob. (1 - q) 

At with prob. \ - (1 - q) 

At with prob. \ - (1 - q) 

r A~i with prob. (1 - q) 



(3.10) 



where / = 



\-{\-q)k 2 

\ i-a-q) ' 



and where k is chosen such that kbs/At matches the (1 - <7)-quantile of the random variable 
y(t+ At)- E t [y(t + At)] for the time-step At. In particular, the distribution for y{t + At) will 
be very close to a Gaussian for small At, and we can compute k from the Gaussian distribution 
function 0() as k = ® _1 (g). For example, if q = 0.995, then k = 2.58 and I = 0.971 (rounded to 
three significant digits). 

Given the quadrinomial discretisation (3.10), we can now proceed as in the previous sections, 
and derive a pde for the price operator n c {t, y): 



n c t + an c v + \b 2 n c vv + 6kb\n c v \ - m c = 0. 



(3.11) 



This pde is exactly the same as (3.4), except for the factor 8k instead of /3 in front of b\n c y \. Of 
course, this should not come as a surprise, since for a small time-step At the (1 - g)-quantile 
of y{t + At) converges to k times the standard deviation b\/At, and hence the cost-of-capital 
pricing operator n c {) should converge to the standard deviation pricing operator 7r s () with (} - 
6k. 

If the payoff f{y[T)) is monotonous in y[T), we can represent the cost-of-capital price n c {t, y) 
in the same way as the standard-deviation price (3.6) with respect to the "risk- adjusted" process 

y 

dy={a{t,y)±8kb{t,y))dt + b{t,y)dW. (3.12) 



3.3 Davis Price 

So far, we have established that the continuous-time limit of the iterated variance price n v {) 
solves the non-linear pricing pde (2.16). The continuous-time limit of the standard deviation 
principle n s Q (and also the Cost-of-Capital principle) solves the linear pricing pde (3.5), which is 
considerably easier to solve. In this section we provide a connection between these two pricing 
principles. We show that the linear standard deviation price can be interpreted as the "small 
perturbation" expansion of non-linear variance price. The core of this idea can be traced back 
to Davis (1997). 

Let us consider the small perturbation expansion in more detail. Suppose we already have 
an existing portfolio of insurance liabilities, where the variance price n v (t, y) has already been 
determined for all relevant t and y. Consider now a small position in an additional insurance 
claim with payoff £g(y(T)) at time T. Let us assume that for small e the total price of the insur- 
ance portfolio can be decomposed into n w {t, y) + en D {t, y), where n D [) denotes the "Davis-price" 
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of the additional claim eg(). The total price n v Q + en D Q should solve the pricing pde (2.16), and 
we find 

in] + £7r°) + ain y + en®) + \b 2 in N yy + en® y ) + 



i r 



2 X e rt 



b 2 [in y ) 2 + 2en v y n® + £ 2 (tt°) 2 ) - r(n v + en D ) = 0. (3.13) 



By definition, the price n v Q solves the pde, and we simplify the expression to 

tt° + an G v + \b 2 T& v + \—^-b 2 hnln^ + ein®) 2 ) - m D = 0. (3.14) 
r y 2 yy 2 Xne rt \ y y y ) 

For small e, we can ignore the non-linear £-term, and we obtain a linear pde for the small per- 
turbation Davis-price 



a + Y-^- t b 2 n v y ^n^ + lb 2 7i^ y + r7T D = 0. (3.15) 



Using the Feynman-Kac formula, we can represent the solution to (3.15) as: 

7r D (?,y) = E? [/(y(T))|y(fl = y] , (3.16) 

where [] denotes the expectation at time t with respect to the "risk- adjusted" process y de- 
fined as 

dy= [a{t,y) + — Y — b 2 {t,y)n y (t,y))dt + b{t,y)dW I) . (3.17) 

Note, that the Davis-price defined only "relative" to existing portfolio price 7r v (). In particular, 
we can interpret the adjustment to the drift term as the risk-aversion y/Xoe rt times the local 
"standard deviation" of the existing price process bit, y)n y {t, y) times the "standard deviation" 
b{t, y) of the insurance process y the drives the additional claim g[yiT)). 

4 BSDE's and g-expectations 

The non-linear pde's we have derived in (2.9), (2.16), (2.19) and (2.29) have been studied in re- 
cent years in the context of backward stochastic differential equations (or BSDE's). Necessary 
conditions for the existence and uniqueness of solutions of BSDE's have been established. For 
an overview of applications of BSDE's in finance, we refer to El Karoui et al. (1997) and Barrieu and El Karoui 
(2009). 

Using "BSDE notation" one can show that the solution to the pde (2.16) can be represented 
by the triplet of processes iy t , Y t , Z t ) satisfying 

' dy t = ait,y t )dt+bit,y t )dW t 
dY t = -git,y t ,Y t ,Z t )dt+Z t dW t (4.1) 

iW(ym), 
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with "generator" g{t, y, Y,Z) = | x J grt Z 2 - rY. The realisation of the stochastic process Y t (de- 
pending on y t ) is then considered to be the solution to (2.16). Another example: to represent the 
solution to (2.29) of the price under the mean value principle for a given function vQ, we would 
use a BSDE with generator g ( t, y, Y, Z) = \ ( v" ( F) I ' V \Y)) Z 2 . 

The mathematical setup of BSDE's allows significant generalisations over the restrictive Marko- 
vian diffusion setting we are using in this paper. For example, one can extend the setup to Levy 
processes, see Nualart and Schoutens (2001). However, we will not explore that avenue in this 
paper. 

Furthermore, solutions to BSDE's are always time-consistent, and can therefore be used to 
define non-linear "g-expectations" (where the prefix "g" refers to the generator g{t,y, Y,Z) in 
(4.1)) and the related notion of "g-martingales", see, for example, Peng (2004). The connection 
between g-expectations and time-consistent risk measures has been explored in Rosazza Gianin 
(2006), Barrieu and El Karoui (2009), Delbaen et al. (2010). 

Link with coherent and convex risk measures, see RosGia. 

5 Good Deal Bound Pricing and Model Ambiguity 

In this section. 

6 Summary and Conclusions 

In this paper we have investigated well-known actuarial premium principles such as the vari- 
ance principle and the standard-deviation principle, and studied their extension into both time- 
consistent and market- consistent directions. The method we used to construct these exten- 
sions was consider one-period valuations, then extend this to a multi-period setting using the 
backward iteration method of Jobert and Rogers (2008) for a given discrete time-step At, and fi- 
nally consider the continuous-time limit for At — >• 0. We showed that the extended variance pre- 
mium principle converges to the non-linear exponential indifference valuation. Furthermore, 
we showed that the extended standard-deviation principle converges to an expectation under 
an equivalent martingale measure. Finally, we showed that the Cost-of-Capital principle, which 
is widely used by the insurance industry, converges to the same limit as the standard-deviation 
principle. 
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